Entanglement and Extreme Spin Squeezing 
for a Fluctuating Number of Indistinguishable Particles 
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We extend the criteria for fc-particle entanglement from the spin squeezing parameter presented 
in [A.S. S0rensen and K. M0lmer, Phys. Rev. Lett. 86, 4431 (2001)] to systems with a fluctating 
number of particles. We also discuss how other spin squeezing inequalities can be generalized to this 
situation. Further, we give an operational meaning to the bounds for cases where the individual 
particles cannot be addressed. As a by-product, this allows us to show that in spin squeezing 
experiments with cold gases the particles are typically distinguishable in practise. Our results 
justify the application of the S0rensen-M0lmer bounds in recent experiments on spin squeezing in 
Bose-Einstein condensates. 
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I. INTRODUCTION 



Spin squeezing [J-Q ^ central concept in quantum 
metrology 0, Q and entanglement detection [g in sys- 
tems with a large number of particles. The most promi- 
nent spin squeezingparameter, defined for N spin-i par- 
ticles or qubits, is [2| 



N{AJi 



(1) 



^ 1 in is a collective spin operator point- 



Here Jn = J2 

ing along the direction n in the Bloch sphere, jn is the 
angular momentum operator for the particle i and _L is 
a direction perpendicular to n. It has been shown that 
a value ^ < 1 implies that the state of the TV particles 
is entangled 0- In addition, it allows for a phase uncer- 
tainty below the shot-noise limit, i.e., A6 < [3,[1,0, 
when used as input of the interferometer implementing 
the unitary transformation e~'^^'^"', where m is a direc- 
tion perpendicular to both n and _L. 

The relation between spin squeezing and entanglement 
has been further extended by S0rensen and M0lmer in 
Ref. [lo| , where bounds on ^ have been derived for a par- 
titioning of the state into groups of at most 1 < fc < iV 
particles. A violation of these bounds implies that there 
is at least one group of more than k particles that is 
fully entangled. Hence the state contains at least (fc-l- 1)- 
particle entanglement or, according to the definition in 
Ref. [l^, an entanglement depth fc + 1. The criteria were 
applied recent experiments on spin squeezing in Bose- 
Einstein condensates (BECs) [ill, However, while 
the criteria were derived for a fixed number of distin- 
guishable atoms, the experiments were performed with 
a fluctuating number of bosons sharing the same trap. 
Hence the criteria have to be generalized to (i) a non- 
fixed number of (ii) indistinguishable particles. 



For the case fc = 1, this has been done in Ref. 
There, it has been shown that 



(Jr.) 



(2) 



is a natural generalization of the spin squeezing parame- 
ter [ij]- In particular, the condition ^ < 1 is sufficient for 
sub shot-noise phase estimation, A9 < —7==, and signals 



entanglement if the input state does not contain coher- 
ences between states with a different number of particles 
[l3j . This justifies the use of the spin squeezing param- 
eter from Eq. (121) in experiments with cold |16l4l8l | and 
ultra-cold lllL [l 2ll atomic gases (for an exhaustive list see 
3). In Ref. [l3l|. it has been argued that, formally, the 
connection between sub shot-noise sensitivity and entan- 
glement holds also for indistinguishable particles. 

In this manuscript, we extend, to the case of a fluc- 
tuating number of particles, the fc-particle entanglement 
criteria of Ref. [l3|. We also show how the generalized 
spin squeezing entanglement criteria of Refs [2l| can 
be extended accordingly. Afterwards, we use additional 
atomic degrees of freedom to extend the fc-particle entan- 
glement criteria to indistinguishable particles. We show 
that in a typical spin squeezing experiment with cold, but 
not ultra-cold atomic gases, the particles can be treated 
as distinguishable effectively. These results apply also 
to other spin squeezing criteria 0, |^, [20l - l30j which are 
generally derived for a fixed number of distinguishable 
particles. 

The article is organized as follows. In Section [Hi we 
discuss the generalization to a nonfixcd number of parti- 
cles. In Section IIIIl we consider the applicability of the 
bounds for indistinguishable particles, discussing explic- 
itly cold atomic ensembles and BECs. The conclusions 
can be found in Section HVl 
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II. SPIN-SQUEEZING BOUNDS FOR A 
FLUCTUATING NUMBER OF PARTICLES 

Let us first recall the definition of (fc + l)-particle en- 
tanglement and how it can be extended to the case of a 
fluctuating number of particles. 

A pure state of N particles is fc-producible [HI, HH if 
it can be written as 



l^fc— prod) 



(3) 



where is a state of 7V„ < fc particles 

(such that J2a=i ~ ^ mixed state is fc- 

producible if it can be written as a mixture p^^p^o^ = 

J:iPi\4%ro,){4%ro,\ with h < k for all I. A state 
that is (fc-l- l)-producible but not fc-producible is referred 
to as (fc-l-l)-particle entangled because it contains full en- 
tanglement of at least one group of fc -f 1 particles (with 
1 < fc < iV) [Hill. The concept of (k-Hl)-particle 
entanglement was referred to as entanglement depth in 
Ref. 

The extension of the above definition to the case of a 
fluctuating number of particles follows Ref. [ll]. In this 
case, we define a quantum state to be fc-producible iff 
it is fc-produciblc in every fixcd-iV subspace. Hence a 
fc-produciblc quantum state without coherences between 
states of different N can be written as 



inc 
r A;— prod 



(4) 



N 



where Pfe^prod is a state of N particles and {Qn} forms 
a probability distribution. In practice, Qat — if TV is 
above some threshold due to energy restrictions in the 
lab. For a general state p which may contain coherences 
between different N, wc introduce the projection 



Iwpljv = QnP 



{N) 



(5) 



where Ijv is the projector to the subspace of N parti- 
cles and p*^^-* is a state on this subspace. We may then 
define a state to be fc-producible in general if p'^^ is fc- 
producible for any N. 

Note that there is an ongoing debate about whether 
or not superpositions between states of different particle 
numbers can actually be created (3^ . It turns out that 
since the angular momentum operator Jn = (BnJii^\ for 
any arbitrary direction n, commutes with the number 
operator N = (BnNIn, such coherences do not have 
any effect for entanglement detection with Jn and its 
moments 12311. 



A. Generalizing the S0rensen-M0lmer criteria to a 
fluctuating TV 

Bounds on ^ have been derived for states of A'' spin-j 
particles among which at most groups of fc particles arc 



entangled [l^. The bounds are computed with the help 
of the function [3a| 



F,iX) 



- min(Ajj^)^ 
J P 



(3„> _ 



X 



(6) 



where the minimization is performed over all states p 
of a spin-j particle which fulfill {jn) / j ~ A" for some 
X <E [0, 1] d^. In Eq. (I6|), j±_ and jn are spin operators 
for the single spin-j particle. It is then shown that for 
fc-producible states, the bound 



(AJi)2 > Nj Fkj 



Nj 



(7) 



holds, where Jn = X^i^iJ"'' ^^'^ analogy for J±, as 
introduced above. Hence if the measured values of ( Jn) 
and (AJj^)^ violate Eq. ([7]), then the state is at least 
(fc-l- l)-particle entangled. For j = i, the state allows for 
a smaller uncertainty in an interferometric protocol than 
any fc-producible state. 

Before generalizing these bounds to states with a non- 
fixed A'^, we remark that a different method, based on the 
Quantum Fisher information, to detect (fc + l)-particle 
entanglement for a state of a fixed number of particles, 
has been recently introduced in Ref. j37| . 
Observation 1. For k-producible states of spin-j parti- 
cles with a fluctuating total number, and with given av- 
erage values (TV) and (Jn); the inequality 



(AJ^)2 > (TV) J Fk, 



( (Jn) \ 

\{N)jJ 



(8) 



holds, irrespectively of whether or not coherences between 
different numbers of particles are present in the state. 

The proof is given in Appendix A. Note that Eq. (I8|) 
reduces to Eq. ([7|) for a fixed number of particles. Also 
for a fixed TV, Observation 1 extends the seminal result 
of Ref. [l^l in two ways. Firstly, in the proof, a step 
is carried out [below Eq. ([26]) ] which was not discussed 
explicitly in the original proof. Further, Observation 1 
does not require ^ to be an integer as in the original 
criterion. In order to apply it for non-fixed TV, simply TV 
has to replaced by (TV), as in the usual spin squeezing 
criterion Il3l. 



B. Generalizing other spin squeezing inequalities 
to a fluctuating N 

We now consider other spin squeezing inequalities for 
entanglement detection d, 0, [20l - [30| . which have been 
derived for a fixed number of particles. Most of them 
can be generalized to the case of a fluctuating number of 
particles by directly using the inequality 

(AJn)^ = E^?^((^''^)')-(EQ^(^^^))' 



N 

> ^Q^(AJW)^ 

N 



N 



(9) 
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which can be derived using the Cauchy-Schwarz inequal- 
ity, for any arbitrary direction n. It also follows from 
the concavity of the variance. Further, note that for any 
operator O = ®'^=qO'^^^ i which commutes with TV, 



{&) = Tr[p6'] = ^Q^Tr[p(^)(0(~))'] 



(10) 



N 



holds for any power I. All angular momentum opera- 
tors Jn are of this form. Therefore, coherences between 
states of different in p do not play any role in entangle- 
ment detection with any kind of spin squeezing criteria, 
as mentioned above. 

As an example, we perform the generalization for the 
complete set of inequalities from Ref. [1^ and for the 
criteria detecting /c-particle entanglement from Ref. [llj . 
All these criteria have been derived for N particles with 
spin i — 2 . 



^ The set of criteria of Ref. [20| is 



( J^) + (J2) + (J2) < 7V(A^ + 2)/4 (11) 
(AJx)' + (A,/y)2 + (AJ,)2 > NI2 (12) 
(Jf ) + (./j') - A^/2 < {N-\){^,hf (13) 
(Af-l)[(AJi)2 + (AJj)2] > ( J2) + iV(A^ - 2)/4, 

(14) 

where i, j, 1 take all possible permutations of x, y, z. This 
set is complete in the sense that it detects all entangled 
states which can be detected based on the knowledge of 
{Jf) and (AJi)2 for i = x,y,z [2|. 

Due to linearity, inequality pT|) . which is valid for all 
quantum states, directly generalizes to 

(ix> + (iy) + (./z><[(^'>+2(TV)]/4. (15) 

Inequality (fT2|) can be generalized using Eq. Q to 



(A4)2 + (AJy)2 + (A J,)2 > (A)/2. 



(16) 



In analogy, the inequalities (jl3|14p can be generalized 
by applying Eq. ([9]) to the variances. The result can be 
written as 



(AJi)" > ((Ar-l)-ijf^ 



(17) 

+ ((iV - 1)^1 Jf) - ((iV - l)-iiV)/2 



(AJO^ + (AJj)^ > ((A-l)-V,-^) 



(18) 



+ ((iV-l)"^A(iV-2))/4. 



Here it is assumed that Qq ~ Q\ = 0. This should not 
pose a problem because the spin squeezing criteria arc 
developed for a large number of particles. A conceptual 
change in the generalized criteria from Eqs (fT7|) and ([18]) 
is that instead of the expectation values ( J^), terms such 
as {{N — appear. This implies that the number 

of particles has to be measured in each shot, which might 
complicate the application in some experiments. In the 
same way, the set of inequalities for A^ spin-j particles 
from Ref. [s^l can be generalized to a non-fixed N. 



Note that alternatively, the criteria could be tested for 
a fixed number of particles A^. In this case, one could col- 
lect separate statistics for each A^. If the number fluctu- 
ates strongly, it would be very difficult to collect enough 
statistics for a given fixed A^, while it is still possible to 
have enough statistics for the generalized criteria. 

We finally remark that the bound 



(AJ.)2 > 



1 



N 



(4) 



N 



1 

4 



(19) 



for fc-producible states from Ref. [2l| can be generalized 
to 

(Air > ^ [{N-'Ji) + {N-\J^)] - 1. (20) 

This bound is optimal for the symmetric twin-Fock states 
with A^/2 particles in each of the two modes of an in- 
terferometer which promises a phase uncertainty close 
to the ultimate Heisenberg limit, A9 — |33|. Re- 
cently, such states have bee n prep ared experimentally 
with ultra-cold atomic gases |39l442l | . Since the number 
of atoms fluctuates in these experiments, Eq. ([20]) could 
be used to bound k, while Eq. ^ from Observation 1 
is generally not useful in this situation since (Jz) ~ 
for these states. However, the same problem concerning 
the indistinguishability of the particles occurs also here. 
This problem will be discussed in the next Section. 



III. SPIN-SQUEEZING BOUNDS FOR 
INDISTINGUISHABLE PARTICLES 

The bounds ([7|) and jS]) presented above have been 
derived for distinguishable particles. This corresponds 
to the usual situation employed in quantum information 
theory with, for instance, trapped ions. In this case the 
particles are assumed to sit at remote locations and oper- 
ations are only performed on the internal degrees of free- 
dom, locally at each trap. The particles can be treated 
as distinguishable, labelled by the trap number, and the 
(anti-)symmetrization can be dropped p^ . 

However, Eq. ([7]) has been recently applied to discuss 
spin squeezing experiments with Bose-Einstein conden- 
sates [ill, mi • In this situation, all the particles (bosons) 
share the same trap state. Their collective internal state 
has to be fully symmetric with respect to the interchange 
of any two particles in first quantization. For indistin- 
guishable bosons, the (symmetric) fully separable states 
have the form \ . The spin squeezing condition ^ < 1, 
see Eqs. ([1]) and ([2|), still holds and signals entanglement 
in the sense that the state of the indistinguishable bosons 
cannot be written as \(j>)^^ . The relation between shot- 
noise limit and separable states holds formally as well 

[H, El HIS- 

In contrast, a symmetric state of A^ particles can be 
either fully separable or fully entangled, but no symmet- 
ric states that are fc-particle entangled as in Eq. ^ exist 
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for 1 < fc < iV [50l - l52| . Hence the classification intro- 
duced above for distinguishable particles is not directly 
applicable to recent experiments with BECs, where the 
individual particles are not addressable. The same prob- 
lem would occur if the criteria for fc-particle entanglement 
proposed in [2lj and generalized in Eq. ([20)) were applied 
to the twin-Fock states produced recently with ultra-cold 
atomic gases [sol - ii^ . 



A. Entanglement and spin squeezing due to 
symmetrization 

First, let us notice that the collective spin operators 
Jn, which appears in the definition of the spin squeez- 
ing parameter ^ [Eqs. ^ and are permutationally 
invariant, i.e., P^J-aP-K = Jn for any of the A^! permuta- 
tions TT of the N particles (represented by Ptt)- Therefore, 
Tr[pJ„] = Tr[ppiJ„], where pm = w\T.7tHpP^ is per- 
mutationally invariant. One may think that, because of 
this property of the collective spin operators, the spin 
squeezing bounds for non-symmetric and corresponding 
symmetrized states should remain the same. However, 
a state of N bosons needs not only to be permutation- 
ally invariant, but symmetric with respect to the inter- 
change of any two particles, i.e., it has to be possible to 
write it as a mixture of symmetric pure states fulfilling 
-P-TT |'0s) = I'/'s) for any permutation tt. This is a much 
stronger requirement [43|. 

Consider, for example, the permutationally invariant 
state of TV = 2 particles gpi = [|0)(0| «) |1)(1| 

|0)(0|]/2. Here |0) and |1) are eigenstates of the Pauli ma- 
trix cTz with eigenvalue 4-1 and -1, respectively. This can 
be rewritten as gn = [|^+)(^"''| + (^Z'" |]/2, where 
= (|01) ± |10))/\/2. Hence gn does not live on 
the symmetric subspace because it has an antisymmetric 
component \ip~). Since the state is separable, ^ > 1 for 
any combination of the directions n and _L. Hence it does 
not allow for sub shot-noise phase estimation. Project- 
ing gpi onto the symmetric subspace leads to gpi — > \ip^). 
This state is known as a twin-Fock state of iV = 2 parti- 
cles [1^ . It is entangled [i^ and allows for sub shot-noise 
phase estimation [38| even though it is not spin-squeezed 
because {ip~^\Jn\il^^) = for any n. In Appendix B, we 
consider an additional example where a separable state 
is transformed into an entangled spin squeezed state by 
symmetrization. 

This shows that symmetrization does not preserve nei- 
ther entanglement nor spin squeezing. In general, sym- 
metrization does not preserve the fe-producibility class 
of a state of N particles. A fc-producible state will gen- 
erally be A^-particle entangled after the symmetrization, 
and the bounds for a given k do not apply anymore. 



B. Generalizing the S0rensen-M0lmer criteria for 
indistinguishable particles 

We assume that the collective spin transformations and 
measurements are performed on two energy levels of each 
atom, which we refer to as the internal degrees of free- 
dom. The extension of the bounds ([7]) and ([5]) to in- 
distinguishable particles is based on the inclusion of the 
atomic external degrees of freedom such as the spatial 
trap states. We thus consider operations of the form 
Ain (8) lox, where Ai^ acts on the internal degrees of free- 
dom and lex is the identity acting on the external degrees 
of freedom. The operator A-^ (8) lox must be permuta- 
tionally invariant because we consider indistinguishable 
particles [i^. As mentioned above, this is the case for 
the collective spin operators Ain = Jn- 

The basic idea is that particles can be distinguished 
here by their external state. Therefore, the state needs 
to be symmetrized only with respect to all particles in 
the same external state, but not with respect to particles 
in different external states. This is true even though the 
operations introduced above do not resolve the external 
states [ssj . 

Let us illustrate this with an example. We consider 
N ~ 2 particles, labeled as 1 and 2, in two different 
external states, labeled as a and b {{a\b) = 0). Following 
Ref. (4^ . a general pure symmetric state can be written 
as 

IV-) = -^{\lpl2}in ® |ai&2)cx + \lp2l)in ® |fcia2)cx), (21) 

where |'!/'i2)in is a general (not necessarily symmetric) 
internal state of the two particles, |'02i)in ~ An|V'i2)in, 
Pin permutes the particles, and \aibj)cx is the external 
(e.g. spatial) wave function — 1,2, i ^ j). The 
mean value of the operator Ai^ (g) lox is 

(^|ii„« loxl^) = + (22) 

where the two terms in the sum are equal since v4.in is 
permutationally invariant: 

(V'2l|ii„|V^2l) = {A2\PLAnPin\i'l2) - (^12 lAn I ^12) • 

(23) 

We dropped the label "in" of |'i/'i2)in for simplicity. 
The above equations show that IV'12) is sufficient to de- 
scribe the state of the two particles. In particular, non- 
symmetric states IV'12) are allowed and the two particles 
can be formally treated as distinguishable. 

The generalization to a system of Nj particles in the 
external level 7 (such that ^7 = ^) can be formu- 
lated as follows. 

Observation 3. The expectation value {Ai^ (E) lex) for 
any permutationally invariant operator Ai^ with respect 
to a fully symmetric state \tJj) with Nj bosons in the 
external state 7 is equal to the expectation value (^in) 
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FIG. 1: Illustration of Observation 3 with particles of din ~ 2 
internal states. The light-gray structure (color online) indi- 
cates entanglement between particles, (a) Fully distinguish- 
able particles (dcx = N). Depicted is a 3-particle entangled 
state of A*' = 8 particles. Entanglement is present between 
the particles with energy level 7 equal to 2,3, and 5 (counting 
upwards from the lowest level with 7 = 1), and between par- 
ticles 7 and 8. (b) Fully indistinguishable particles (dox ~ 1). 
Depicted is a state which does not factorize. (c) Mixed sit- 
uation (1 < dcx < N). Depicted is a 4-particle entangled 
state of = 11 particles in dex = 5 external levels. There 
are two groups of 4 fully entangled particles: in the external 
level 7 = 2, where the 4 particles are indistinguishable and in 
a non-factorisable state; and in levels 4 and 5. In the latter 
case, the group of particles in level 7 = 4 is distinguishable 
from the particle in level 7 = 5. The two particles in the 
lowest level 7 = 1 are in a symmetric separable state of the 
form \(p) (g) 10). 



for particles that share the same external level 7. In this 
case, the N-y particles may be only found in a fully entan- 
gled or fully separable state. However, entanglement is 
also possible between particles occupying different levels. 

We can now extend the S0rensen-M0lmer bounds. A 
state can be considered as (effectively) fc-producible if 

|^)i„ = ®*iilV'a). (24) 

where \ipa) is state of Na < k particles (X^aLi = N) 
for all a. The particles in the state \tpa) can occupy 
a single external state 7 (in which case Na = and 
IV'q) = I "07) is symmetric) or different external states 
J G la (in which case Na ~ J2-yeia ^^'^ 1^"'* 
not necessarily symmetric). As an example, the state 
schematically shown in Fig. [T](c) is 4-particle entangled. 

With this notion, the S0rensen-M0lmer criteria can be 
applied in systems of indistinguishable particles as fol- 
lows. 

Observation 4. (i) For particles of spin j, if the spin 
squeezing parameter violates Eq. ([^ for a given k, then 
the input state cannot be written as a mixture of effec- 
tively k-producible states of Eq. {24^ . (ii) For particles of 
spin i, if the spin squeezing parameter violates Eq. ^ 
for a given k, then the input state allows for a smaller 
phase uncertainty than the smallest one achievable with 
a mixture of effectively k-producible states of Eq. \24^ . 

In both cases, effective {k + l)-particle entanglement 
is proven by a violation of the criteria. These notions 
directly generalize to systems of a fluctuating number of 
particles as in Section HH 

C. Cold atoms 



computed with respect to the corresponding internal state 
\il^)i-a, which is symmetrized only with respect to the par- 
ticles sharing the same state 7, for all 7. 

This Observation is formulated more precisely in Ap- 
pendix C below, where also the relationship of \tp) and 
|'(/')in is explained in detail. 

Figure [T] illustrates several examples of particles in 
din = 2 internal and dex external modes. The usual sit- 
uation employed in quantum information theory, where 
all particles are distinguishable (dox = N), is shown in 
Fig. [T](a). For what concerns our discussion, this is for- 
mally equivalent to an array of separated wells, as in ions 
traps. The opposite situation of all particles occupying 
the same level (dox = 1), is shown in Fig. [TJb). For 
indistinguishable particles, only two possibilities are al- 
lowed in this case: either all particles are in a separable 
{i.e., product \(j>)'^^) state, or all p articles are entangled, 
due to the symmetrization |50l452| . As mentioned above, 
the fc-particle entanglement criterion discussed in Sec. |TT] 
does not apply in this case. The interesting intermediate 
situation is shown in Fig. [T](c). In this case several par- 
ticles may occupy the same external state. As noticed 
in Observation 3, the symmetrization is necessary only 



Observation 3 is also useful in the context of entan- 
glement detection with generalized spin squeezing in- 
equalities (SSI) d, Q in cold atomic clouds. Usually 
the atomic ensembles are not ultracold, and can be as- 
sumed to be in a thermal state externally. We estimate 
the population of the trap levels using the statistics of 
an ideal Bose gas taking the parameters from a, typical 
experiment with a cigar-shaped configuration (54 . [ssl ]: 
ujz = 2, uj± — 1000 (trap frequencies in units of 27rs~^), 
T = 30^K, = 5 • 10^ The chemical potential fi is 
defined implicitly by the relation N = '^/^{rik), where 



{Tlk) = 



-fin(Ek-ti.) 



(25) 



is the average population of the level k, /3 = l/ksT with 
the the Boltzmann constant fcs, and is the energy of 
level k. We approximate this by the energy levels of a 3d 
harmonic oscillator with the given trap frequencies. 

The largest average population is obtained for the 
ground state with energy Eq = h{ujz/2 + uj±). We obtain 
l3fi sa -12.41, which leads to (no) « 4.1 • 10"^. Further, 
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the probability of having particles in the level k, given 
by the expression in square brackets in Eq. p5|) . decreases 
exponentially with rifc . The ratio Pn^+i/puk = e~^^^''~^'> 
is the largest for k = 0, where it is equal to 4.1 • 10~^. 
Therefore, it can be assumed that at most one particle 
occupies each level. Since only internal quantities are 
used in the generalized SSI, one can therefore treat the 
particles as distinguishable by using Obs. 3. This situa- 
tion corresponds to the one depicted in Fig.[lja). 



making use of additional degrees of freedom of the atoms. 
The spin squeezing bounds of Ref. [l^ can then be inter- 
preted as conditions of such effective fc-particle entangle- 
ment. Our results make it possible to apply the bounds of 
Ref. (lo| in spin squeezing experiments with cold atoms 
and Bose-Einstein condensates. 
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D. Bose-Einstein condensates 

In Bose-Einstein condensates particles share the same 
external state, as schematically illustrated in Fig. [ijb). 
In this case, it is not possible to directly apply the bounds 
(I?!) and ([5]): particle entanglement is either absent or 
maximal, due to symmetrization (sol - fs^ . 

However, we recall that Eqs. ([7]) and ([8]) are sufficient 
conditions for entanglement and, as mentioned in the in- 
troduction, for spin j — \ particles, they are strongly 
related to the usefulness of the entangled state for pa- 
rameter estimation. Therefore, a measurement of the 
spin squeezing bounds in BECs might still be consistent 
with k < N , even if, by some other means, it is possible 
to show that all particles share the same external state, 
as in Ref. [llj ■ The outcome k < N should be interpreted 
either as due to noise or by saying that all particles are 
entangled but the state is only partially useful for pa- 
rameter estimation. It is known, indeed, that there are 
symmetric states which are fully TV-particle entangled, 
but that are not spin squeezed and do not allow for sub 
shot-noise phase estimation [4^ . 

Finally, one might think that, by making the BEG 
cloud very dilute, it is possible to effectively distinguish 
the particles and thus use the spin squeezing bounds (III) 
and ([5]). In Appendix D we show that this is not the 
case. 



IV. CONCLUSIONS 

The spin squeezing criteria introduced by S0rensen and 
M0lmer for N distinguishable particles in Ref. [l^l are a 
powerful and experimentally feasible method to detect k- 
particle entanglement, also referred to as entanglement 
depth k. However, most of the spin squeezing experi- 
ments are performed with a fluctuating number of parti- 
cles and, as in the case of BEG, these particles are indis- 
tinguishable. To fill this gap between theory and experi- 
ment, we have extended, in the first part of this article, 
the S0rensen and M0lmer criteria to systems with a fluc- 
tuating number of particles. We have also shown how 
other spin squeezing inequalities [lO, |2l[ can be general- 
ized to this situation. In the second part of the paper, we 
discussed the conceptual problems that occur when the 
individual particles are indistinguishable. In this case, 
effective fc-particle entanglement can be defined only by 
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APPENDIX 

A. Proof of Observation 1 

The proof follows the lines of the proof of Eq. ^ for 
fixed N [l^l using methods developed for non-fixed N in 
Ref. [l3| . We want to compute a lower bound on the vari- 
ance of J]_ for all fc-producible states. Since the variance 
is concave in the state, its minimum value is reached by 



pure states of the form IV'fc-prod) =YjN V^lV'i^prod) 
IstI . where \/Qn are real numbers with ^j^Qn = 1 
i„;,(^a)\ fc-producible state 

Using Eq. ^ we can 

(^^)^2 



and l^lZlo,) = 1^^" 
of N particles [cf. Eq. ©]. 

write (AJj_)2 > J2n QNi^J^ ')^- 1^ addition, we 
note that, due to the product structure of the states 
l^fc-prod) and since is the sum of operators act- 

ing on fixed Na subspaces, = J2t!=i J±^''\ have 
(Ajf'V = Ea=i (Aj|^°y . Therefore, the variance 
of J± for fc'-produciblc states is bounded by 



Mis 



(26) 



N 



where the operator on the right hand side acts on the Na 

particles in the state iV'o^"'')- Note that we did not attach 
an index N to Na in order to simplify the notation. 

Now we have to find the minimal bound for the vari- 
ances (Ajf'°')2 for every N, given the mean value 

(J^"'). If we consider 7V„ spin-j particles, the total spin 
ja can range from (if Na is even) or 1/2 (if Na is odd) 



the small- 



up to Naj- Wc show that for any (j, 
est bound is reached by choosing the largest total spin 
possible by using that (Aj|^°^)2 > jaFj^{{jL^"^) /ja) 
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for states with a fixed spin ja, cf. Eq. ([6]). The in- 
gredients needed for this, which have been proven in 
Ref. are: (i) the functions Fj^ are convex, i.e., 

Fj^ [aX+bY) < aFj^ {X)+bFj^ {Y) for aU j„ and a,b>0 
with a + b = 1, (ii) Fj^ (0) = for aU ja, and (iii) that 
Fj^iX) < Fj,^{X) if "ia > fa- By using the inequal- 
ity (i) with a = Y = and b = 1 — a, and the 

property (ii), we have FjJf^X) < f^Fj^iX) / j^. Mul- 
tiplying by ja both terms and using (iii), we arrive at 
jaF,^[^X) < j'aF,,jX) if fa < Ja- Finally, taking 



X 



, we have 



B. Example: symmetrization creates spin squeezing 



Let us consider the state 



IV'a) = V^|ii) + VT^|oi) = (V^|i)+VT~~^|o))®|i) 

(30) 

which is clearly separable. Therefore, it has a spin 
squeezing parameter > 1 for any a and any combina- 
tion of the directions n and _L, i.e., it is not spin squeezed. 
In particular, for the directions n = z and _L = x it is 
given by [cf. Eq. (P)] 



1 



2 > 1. 



(31) 



jaFj^ 



Ja 



< JaFf^ 



fa 



(27) 



if ja 5^ j'a- Let us now consider a superposition lip) = 
Cj^\ipj^) + Cj^lV'j^) of states with a different fixed spin 
ja > j'a ■ Since the spin operator J± does not couple the 
states of different total spin ja and ja , its variance with 
respect to lip) is equal to the variance with respect to the 
mixture p = \cjj'^\ ipj^ ){ipjj + \ Cj^^ P | V-j J {'fpj'^ \ ■ Using 
the concavity of the variance, we obtain that 

> |Cjj'f;.((in)|^^)/ja) + \Cj'fF,,J{Jn)^^^/fa) 

> Fj^{{Jn)l^)/ja), 



where we have used Eq. (f27|) and |c 



1. 



Taking the maximum value of ja (i.e., ja = Naj) we 
arrive at 



The corresponding state which could be realized in the 
scenario of indistinguishable bosons is obtained by sym- 
metrizing (and normalizing) the state \ipa)-, leading to 



|oi) + |io) 

V2 - 



(32) 



where /? = iq^- For this state, the spin squeezing param- 
eter for the same particular directions n = z and ± = x 
is given by [cf. Eq. ([T])] 



4-4. 



(33) 



This is smaller than the critical value 1 for a larger range 
of parameters. The minimum is reached at = 7/8 
for /3 = 4/5, which corresponds to a = 2/3. Therefore, 
symmetrization does not preserve neither entanglement 
nor spin squeezing. In order to illustrate the results, the 
two curves are plotted in Fig. [2l 



AJ 



> NajF, 



Naj 



> NajFkj 



Naj 



(28) 

where the second inequality is due to (iii) and Na < k for 
fc-producible states. Since the function Fj^ (^IJ^ convex 
in X, we can now apply Jensen's inequality [59[ to the 
last term in Eq. p8| . We obtain that 



Mjv 



a=l 



(JV)\ 



Nj 



(29) 



J^ai"^^"') f°^' N. Finally, by com- 



where (J, 

bining Eqs. p6|) and ((29|) , and using again Jensen's in- 
equality [59|, we have 



(JA' 



> {N)jFkJ^], 



where (N) = Eiv QnN and (J„) = Qjv(J^'^'). This 
proves Obs. 1 [cf. Eq. dS])]. □ 



C. Proof of Observation 3 

Let us first introduce the formalism used in the proof. 
Recall that we are considering here N particles with c?in 
(dcx) internal (external) degrees of freedom, labeled by 
i = l,...,din (7 = l,...,dox)- Referring to Fig. (P), we 
can think of the external states as the energy levels of 
a spatial trap. Each level 7 contains iV^ particles in the 
state |7)ex- These particles can be in different internal 
states |i)in. In general, we have Ni^^ particles in the state 
\i,j) = \i)in <8) |7)cx (to simplify the notation, we remove 
here the tensor product sign and pendices "in" and "ex" ) . 
We also introduce a vector Nex = {Ni,N2, ...,Nd^^) giv- 
ing the occupation numbers of each external state and a 



vector 



2,7, 



Ndi^,^) with occupations of 



the internal states for a fixed external level 7. Here, the 

relations J2j=i = N and X]f=i -^i,7 = -^7 hold. 

Let us consider a specific example for c?in = dcx = 2 
and iV = 3. Choosing TVi.i = 1, N2,i = 0, Ni-2 = 2 
and N2,2 = 0, wc obtain Nox = (1,2)^ Ni = (1,0) and 
N2 = (2,0). The (non-symmetric) state is 



^7=1 



^LiK,!)^"^"^ = |1,1)®|1,2)®|1,2). (34) 
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FIG. 2: Plot of the spin squeezing parameters [dashed 



hne, Eq. ((3T])] and [sohd line, Eq. dSS}]. While never 
goes below the value 1 [dotted line], is below 1 for a large 
interval. 



Finally, {N^}^ is the complete set of occupation numbers 
N-y for all the 7 levels. The corresponding symmetrized 
states with occupation numbers is given by 



\D 



{N^}\ ^ 



1 



'i=i\hi) 



(35) 



where i-^ is a representation of the permutation tt, and 
the sum runs through all distinct permutations, the num- 
ber of which is A/" = ({n' }) = 1 . The states 



\D- 



N„x / 



form a basis which is the analogous to a Fock 
state basis in second quantization. 

The non-symmetric state from the example above [cf. 
Eq. ([34])] becomes 



\D- 



{N-,}^ 



+ |1,2) ® 



® |1,2) ® |1,2) -I- 

|1,2) + |1,2) 



track of the external level 7 the particle is in. This will 
be important below. 

With these definitions, we can reformulate Observation 
3 in technical terms. 

Observation 3. For any permutationally invariant op- 
erator acting on the internal degree of freedom, and 



for a symmetric state j^'g"") = X]{n }'^{n^}I-C' 
with a fixed occupation vector Ncx, 



Ne 



Ne 



holds, where 



Ne 



(38) 



I^{N^}C{N^}[ ®7 1^^;)]' 



rN-, 



|/^^) is a symmetric internal state as defined in Eq. (5^ 



Proof. By inserting the definitions of 1^'^°'') and 



I* 



Ne 



it is easy to see that Eq. ([38]) holds if 



Nex / 



In:\] a 



V,)] 

■y 

(39) 

is true for all {N^} and {Ny} with the same Ncx. We 
will show now that this is the case. We insert into the left 

Nex ' 



hand side of Eq. ([5^ the definition of the states \D- 
[cf. Eq. ([35]) ]. which leads to 



E 

tt.tt' 



t . 
A, 



l« ,7 



(40) 

As before, we consider the sum of distinct permutations 
only. Due to the identity lox on the external states the 
terms in the sum will vanish unless the Nj particles in 
level 7 are on the same positions in the permutations 
on both sides of Ai^ (g) lex, since (i,7|Ain (g) loxK',7') = 
{i\Ain\i')S^ 
as 



Therefore, we can rewrite expression ([401 



i:[ 



p p- 



Af 



P P-/ 



IN' 



(41) 



We will use the label D in general for Fock states with 
a fixed occupation in internal and external levels in first 
quantization. In particular, we employ symmetric states 
with particles in the single external level 7, 



\D- 



\I 



Here the permutations tt and tt' permute particles stem- 
|1, 2) (g) |1, 1)) . ming from the same external state 7, and tt permutes 
particles with a different 7. Note that for simplicity we 
use the same operators to represent a permutation tt 
of the N particles even though now the state space of 
each particle is reduced to the internal states. In order 
to clarify the notation we employed, we note that for the 
example considered in Eq. (IM)) . the reduced state would 
be 



where 



(36) 



(37) 



ID 



ID 



ID 



= J = j-j I is the number of distinct per- 
mutations TT. We attach the label 7 to \i) in order to keep 



Since is permutationally invariant, we have that 
P^Ain (g lox^TT = ^in ^ lex- Hcuce in the sum over 
TT each term contributes equally, and the sum can be re- 
placed by the number of distinct permutations Af^x = 
in expression ([51]). 



N 

VNe 



JV! 
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FIG. 3: Initially all particles are in the lowest energy level of 
the trap. The cloud is then released. If it is diluted enough 
along the horizontal direction, it is likely that at most one 
particle falls into each of the boxes, which represent the single 
detectors. 



We observe that because the permutations tt only per- 
mute particles with the same 7, one can rewrite = 
^■y{^7T ^'^i) '^^ permutations tt^ which permute parti- 
cles in the level 7. This leads to 



cf. Eq. (1371). We arrive at 

/aA2 n tv^'N', , N , . , N', T 



One can directly check that the prefaetor is equal to 1. 
Therefore, condition Eq. (l39l) is fulfilled. q 



D. Dilute cloud argument 

One may think that, after preparing the Bose-Einstein 
condensate atoms in the ground state of a confining trap, 
it is possible to apply the spin squeezing inequalities 
and (HI) by simply releasing the trap, letting the cloud 
expand and fall onto a grid of small detectors capable 
of measuring the internal state of a single atom. If the 
cloud is dilute enough, it is very likely that at most a 
single atom enters each detector, thus making the atoms 
distinguishable. This would make it possible to apply the 
S0rensen-M0lmer bounds. The situation is illustrated in 
Fig. 121 We show here that this argument, which is often 
encountered in discussions, does resolve the problem. 

Let us assume that, before releasing the atoms from 
the trap, their state is of the form 



1^) 



|^s)in ® |0) 



cx ' 



(42) 



where |V's)in is a symmetric total internal state and each 
atom is in the ground state of the trap |0)ex- Here "ex" 
( "in" ) indicates the external (internal) degree of freedom 
as in Section IIIIBI We assume here that all atoms share 
the same spatial wave function, which can thus be faetor- 
ized. If interactions can be neglected during free fall, then 



only the spatial state of each atom changes, leaving the 
internal state symmetric. By waiting long enough, the 
single-particle spacial wave function becomes so spread 
that the probability to detect two atoms at the same 
spatial detector is negligible. 

Let us assume for simplicity that the atoms are trapped 
and detected state-insensitively first, such that at most 
one atom is detected in each site. A problem is that in 
each shot, different sites will be occupied. This might still 
be considered as a minor problem. In a one-dimensional 
trap, for instance, it could be resolved by identifying par- 
ticle "1" with the leftmost trap, particle "2" with the 
particle right from particle "1" , and so on. Alternatively, 
one could postselect on events where always the same 
sites are occupied. 

In general, the position measurement makes the state 
effectively distinguishable. Let us illustrate the situation 
with an example for N = 2 particles, labeled as 1 and 2, 
in two different sites labeled by a and b. As in Section 
IIIIBI we consider a general pure symmetric state 



1 



V2 



(1^1 



), (43) 



with the same definitions as in Eq. (l2Tl) . An operator Ma 
acting on the internal state of the particle in site a can 
be written as 

Ma = (il «) ll2)in «) (m? ® a2)ox 
+ (li®i2)in® (ll®m2)ox, 

where m°|a) = \a) and rh°'\b) ~ since we measure lo- 
cally at site a. Ma has to be permutationally invariant 
with respect to the interchange of the particle labels since 
the particles are indistinguishable [4^. The expectation 
value with respect to the state (|43l) is 

ii^lMaH) = l[{lpl2\Ai ® 12|V'12> + (V'2l|ll '»i2|V'2l)], 

(44) 



where the two terms are equal since 

(^/>2l|ll0i2|V'2l) = (^-121^1(11 «'^)A„|^Al2) 
= (^/'12|il «) l2|V'12). 

We dropped the label "in" of |V'i2)in for simphcity. An 
analogous result is obtained when considering an opera- 
tor acting on the internal state of the particle on site b. 
Since only such operators are measured in the usual sce- 
nario, we can identify particle 1 with site a and particle 2 
with site 6, and |'0i2) is sufficient to describe the state of 
the two particles. This is a state of two distinguishable 
particles. 

However, since the measurement acts only on the ex- 
ternal degrees of freedom, the product structure be- 
tween the internal and the external degrees of freedom 
in Eq. (H^l) is preserved. It is then evident that the inter- 
nal state remains fully symmetric even after the position 
measurement, since this affects the external state only. 
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Hence making the state distinguishable effectively after 
the state transformation in the measurement does not 
make it possible to leave the restricted class of symmet- 
ric states: in the interferometric situation we considered 
above, when only the sites a and b are occupied, then 
we arrive at the state of Eq. but with a symmet- 



ric internal state |'0i2)- Therefore, the effective state of 
the distinguishable particles is symmetric. This example, 
which can be directly generalized to N particles, shows 
that simply making the cloud dilute does not make it 
possible to apply the spin squeezing bounds discussed in 
this paper. 
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